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This  paper  is  concerned  with  waves  in  a  composite  elastic  bar, 
the  left  half  of  which  is  composed  of  a  linearly  elastic  material, 
while  the  nonlinearly  elastic  material  of  the  right  half  can 
undergo  a  phase  transition.  We  assume  that  a  wave  in  the  left 
portion  of  the  bar  is  incident  upon  the  interface  between  the  two 
materials,  and  we  investigate  the  question  of  whether  the  phase 
transition  can  be  exploit^  to  augment  or  diminish  the  strength 
of  the  reflected  or  transmitted  wave. 
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1.  Introduction.  We  consider  a  composite  tensile  specimen  consisting  of  two  dissimilar 
elastic  bars  joined  end-to-end.  One  of  the  two  elastic  materials  ("material  2")  is  capable  of 
undergoing  a  stress-induced  transition  to  a  second  phase,  while  the  other  ("material  1")  is  not. 
When  an  incident  wave  traveling  in  the  single-phase  material  strikes  the  bimaterial  interface,  it 
may  nucleate  a  phase  transition  in  material  2.  The  reflection  and  transmission  characteristics  of 
the  product  phase  in  such  a  transition  will  in  general  differ  from  those  of  the  parent  phase.  We 
designate  as  "material  3"  the  single-phase  elastic  material  whose  properties  are  identical  with 
those  of  the  parent  phase  of  material  2.  Our  interest  lies  in  comparing  the  strengths,  relative  to 
the  incident  wave,  of  the  reflected  and  transmitted  waves  generated  in  the  material  1/material  2 


composite  bar  with  the  relative  strengths  of  those  that  would  occur  in  a  bar  composed  of  material 
1  and  material  3.  In  particular,  we  investigate  the  extent  to  which  the  phase-transforming  capa¬ 
bility  of  material  2  can  be  exploited  to  control  the  strength  of  the  reflected  or  transmitted  wave. 


The  composite  bar  is  viewed  as  consisting  of  two  perfectly  bonded,  semi-infinite 
one-dimensional  elastic  continua.  Material  1  is  taken  to  be  linearly  elastic,  while  material  2  is 
treated  on  the  basis  of  a  nonlinear  continuum  model  of  the  kind  currently  receiving  much 
attention  in  discussions  of  the  macroscopic  aspects  of  phase  transformations  of  martensitic  type; 
see,  for  example  the  references  cited  in  [1].  This  model  is  based  on  a  "two-well"  potential 
energy  for  material  2,  and  it  includes  a  nucleation  criterion  and  a  kinetic  relation  governing  the 
initiation  and  evolution,  respectively,  of  the  phase  transition.  Thermal  effects  are  neglected  here, 
so  that  the  model  is  a  purely  mechanical  one.  It  has  been  applied  to  quasi-static  phase  transform¬ 
ations  in  [2],  where  it  was  shown  to  lead  to  results  in  qualitative  agreement  with  some  experi¬ 
mental  observations  on  materials  of  the  shape-memory  type.  The  dynamics  of  the  model  have 
been  studied  in  [3]. 

There  are  many  recent  papers  devoted  to  the  general  issue  of  the  continuum-mechanical 
modeling  of  the  macroscopic  effects  of  phase  transformations;  examples  may  be  found  in  the 
references  cited  in  [1].  The  only  work  of  which  we  are  aware  that  is  related  to  the  specific 
problem  under  discussion  here  is  that  of  Pence  [4,  5].  In  [4],  Pence  studies  the  reflection  and 
transmission  of  an  acoustic  shear  wave  from  an  initially  stationary  phase  boundary  in  an  elastic 
solid.  The  analysis  in  [5]  is  concerned  with  the  structure  of  the  fields  in  two  elastic  bars,  one  of 
which  (the  impactor)  is  composed  of  a  single-phase  material,  while  the  other  (the  target)  is  made 
of  a  material  capable  of  sustaining  a  phase  transformation.  Neither  of  these  papers  addresses  the 
issues  of  primary  interest  here. 

We  describe  the  basic  model  to  be  used  in  the  following  section.  Section  3  contains  the 
formulation  of  the  underlying  wave  propagation  problem,  and  Section  4  is  devoted  to  its 
solution.  In  Section  5,  we  show  that  the  ratio  of  the  relative  strength  of  the  reflected  wave  in  the 
presence  of  the  phase  transition  to  its  relative  strength  in  the  absence  of  the  transition  is  governed 
by  two  material  parameters:  one  is  the  ratio  of  the  mechanical  impedances  of  material  1  and  the 


-2- 


parent  phase  of  material  2,  while  the  other  parameter  is  inherently  related  to  the  phase-transition 
properties  of  material  2.  In  terms  of  these  two  parameters,  we  derive  in  Section  5  conditions 
under  which  (i)  the  reflectivity  is  always  increased  by  the  phase  transition,  independently  of  the 
kinetic  relation  and  the  nucleation  criterion,  (ii)  the  reflectivity  is  always  decreased,  and  (iii)  the 
effect  of  the  transition  on  reflectivity  depends  on  the  details  of  kinetics  and  nucleation.  Some 
analogous  results  for  the  transmitted  wave  are  stated  without  proof  in  Section  6. 

2.  The  model.  In  a  reference  configuration,  the  composite  bar  occupies  the  entire  x-axis, 
with  material  1  in  x  <  0,  material  2  in  x  >  0;  the  referential  cross-sectional  area  is  A.  We  treat 
longitudinal  motions  in  which  a  particle  at  x  in  the  reference  state  is  carried  to  the  point 
x  +  u(x,  t)  at  time  t.  The  displacement  u  is  to  be  continuous  with  piecewise  continuous  first  and 
second  derivatives.  The  strain  and  particle  velocity  are  defined  by  7(x,  t)  =  u^(x,  t)  and 

v(x,  t)  =  u^(x,  t) ,  respectively,  where  the  derivatives  exist;  in  order  to  assure  that  the  mapping 
X  ->  X  +  u  is  one-to-one,  we  require  that  y  >  -1  everywhere.  At  points  in  the  x,t-plane  where  the 
fields  are  smooth,  balance  of  linear  momentum  in  the  absence  of  body  force  requires  that 


=  p(x)  v^ , 


where  a(x,  t)  is  the  stress  in  the  bar,  and 


X  <  0, 

X  >  0; 


(2.1) 


(2.2) 


the  constants  Pj  and  p^  are  the  mass  densities  of  materials  1  and  2  in  the  reference  state, 
respectively.  From  the  definitions  of  y  and  v,  one  also  has 


-3- 


(2.3) 


V 

X 


=  0 


where  the  fields  are  smooth. 


The  stress-strain  relation  for  the  composite  material  1/material  2  bar  is 


a  =  a(Y,x)5^ 


a,(V) . 


x<0,  Y>-1, 
x>0,  Y>-1. 


(2.4) 


where  Oj  and  are  the  respective  stress  response  functions  of  materials  1  and  2  .  For  material 
1,  we  assume  that 


Ol(Y)  =  4ir.  -1<Y<~,  (2.5) 


where  is  Young’s  modulus.  The  stress  response  function  for  material  2  is  taken  to  have  the 
following  "trilinear"  form; 


02(Y)  =  j 


■  ‘"m  ■  ■  TmWt™  ■  y„).  y„  <  y  s  r„ . 


'M' 


'M 


'm 


(2.6) 


p,(Y-Yo).  Y^Y„- 


m 


The  graph  of  O^iy)  is  shown  schematically  in  Figure  1;  the  significance  of  the  constants  Yq- 
JxAy  Y  >  and  a  may  be  read  from  the  graph.  Each  of  the  three  branches  of  the  stress-strain 

M  m  M  m  or 

curve  in  the  figure  is  identified  with  a  phase  of  material  2:  the  two  rising  branches  are  the 
low-strain  and  high-strain  phases,  and  the  declining  branch  is  the  unstable  phase.  For  simplicity, 
we  have  assumed  that  the  Young’s  modulus  p^  is  the  same  in  both  the  low-strain  and  high-strain 
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phases  of  material  2,  but  this  is  by  no  means  necessary.  The  transformation  strain  Yq  is  the 

distance  between  the  low-strain  and  high-strain  branches  of  the  stress-strain  curve  at  any  given 
value  of  stress;  because  a.,  >  a  ,onehasY/.>Y  -Y»*-IfY«<Y  ,  the  stress  o  at  the  local 
minimum  is  positive,  as  in  the  figure.  If  >  y^,  this  minimum  stress  is  non-positive,  and  the 
right  half  of  the  bar  possesses  more  than  one  unstressed  configuration.  The  stress  shown  in 
Figure  1  is  such  that  the  shaded  areas  are  equal;  is  called  the  Maxwell  stress.  It  is  given  by 

for  simplicity,  we  assume  that  >  0. 

Equations  (2.  l)-(2.6)  yield  a  system  of  two  differential  equations  for  y  and  v  for  x  <  0 
corresponding  to  material  1,  and  a  second  pair  of  differential  equations  for  y  and  v  in  the  interval 
X  >  0  for  the  part  of  the  bar  made  of  material  2. 

Suppose  that  there  is  a  discontinuity  in  strain  or  particle  velocity  along  the  curve  x  =  s(t) 
in  the  x,t-plane.  Balance  of  linear  momentum  and  the  assumed  smoothness  of  the  motion  require 
that  the  following  jump  conditions  hold: 

[[CT]]=-s[[pv]],  (2.8) 

s[[y]]  =  -[[v]],  (2.9) 

where  for  any  g(x,  t),  we  vmte  [[g]]  =  [[g(x,t)]]  s  g(s(t)-H,  t)  -  g(s(t)-,  t). 

The  strain  energy  per  unit  reference  volume  for  the  composite  bar  is 
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W(Y,x)  = 


(2.10) 


Wj(Y),  x<0,  Y>-1. 

W^Cy),  x>0,  y>-1. 

where  the  separate  energy  densities  for  materials  1  and  2  are  given  by 

Y 

W^(Y)  =  \(Y')dY',  Y>-1,  k  =  l,2.  (2.11) 

The  potential  energies  for  the  two  materials  are  defined  by  Gj^(y,  o)  =  Wj^(y)  -  cry .  For  every 

fixed  G,  Gj(y,  o)  as  a  function  of  y  has  a  single  minimum  at  y  =  tj/|i.  In  contrast,  the  potential  G^ 

for  material  2  is  such  that,  at  a  fixed  value  of  o  between  g  and  G. .,  G^(-,  g)  has  two  local 

m  M2 

minima  separated  by  a  local  maximum.  When  G  is  outside  this  range,  G2(-,  g)  has  only  one 
minimum.  Thus  when  the  stress  is  such  that  material  2  can  exist  in  either  the  low-  or  the 
high-strain  phase,  is  a  "two-well"  potential  typical  of  two-phase  materials. 

Suppose  that  in  a  piece  of  the  bar  corresponding  to  Xj  <  x  <  x^  ,  the  strain  and  particle 
velocity  are  discontinuous  at  x  =  s(t)  but  are  otherwise  smooth.  Writing 

E(t)  =  J  [w(Y(x,  t),  x) -I- %(x)  v^  (x,  t)]  A  dx  (2.12) 

’‘i 

for  the  total  energy  in  this  portion  of  the  bar  at  time  t,  one  can  show  by  a  direct  calculation  that 

g(X2,  t)Av(x2,t)  -  G(Xj,  t)Av(Xj,  t)  -  E(t)  =  f(t)A  s(t),  (2.13) 

where 
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f(t)  =  [[W(Y(x,t).  X)]]  -  <a(x,t)>  [[Y(x,  t)]], 


(2.14) 


and  <a(x,  t)>  =  (l/2){CT(s(t)+,t)  +  a(s(t)-,t)}.  We  call  f(t)  the  driving  traction  acting  on  the 
discontinuity  at  x  =  s(t).  In  view  of  (2.13),  one  may  think  of  f(t)A  s(t)  as  the  instantaneous 
dissipation  rate  associated  with  the  discontinuity;  if  either  f(t)  =  0  or  the  discontinuity  is 
stationary  (in  the  Lagrangian  sense)  so  that  s(t)  =  0,  this  dissipation  rate  vanishes.  At  any 
discontinuity,  it  iS  required  that  the  dissipation  rate  be  non-negative: 

f(t)k(t)>0.  (2.15) 

One  can  show  that,  if  the  motion  of  the  bar  is  viewed  as  occurring  isothermally,  (2.15)  is  a 
consequence  of  the  second  law  of  thermodynamics;  we  shall  therefore  refer  to  (2.15)  as  the 
entropy  inequality. 

In  the  problem  to  be  treated  here,  we  shall  be  concerned  with  three  different  types  of 
strain  discontinuity.  The  first  of  these  occurs  at  the  bimaterial  interface  at  x  =  0;  because  this 
discontinuity  is  stationary,  (2.15)  is  trivially  satisfied.  The  second  type  of  strain  discontinuity  to 
be  encountered  involves  a  strain  jump  either  in  material  1  or  between  strains  7(s(t)±,t),  both  of 
which  belong  to  the  same  phase  in  material  2.  Such  a  discontinuity  is  a  sound  wave.  Because 
the  stress-strain  relation  is  linear  between  the  strains  7(s(t)-,  t)  and  Y(s(t)-i-,  t)  in  either  of  these 
circumstances,  the  definition  (2.14)  yields  f(t)  =  0,  so  that  (2.15)  is  automatically  satisfied  at  a 
sound  wave  as  well.  Finally,  we  shall  need  to  deal  with  strain  jumps  in  material  2  for  which 
7(s(t)-,  t)  is  in  the  high-strain  phase  while  y(s(t)+,  t)  is  in  the  low-strain  phase.  Such  a 
discontinuity  is  an  example  of  a  phase  boundary,  for  the  problem  to  be  considered  here,  phase 
boundaries  in  material  2  will  always  move  to  the  right,  so  that  s(t)  >  0.  It  then  follows  from 
(2.15)  that  f(t)  must  be  non-negative  at  the  phase  boundaries  arising  in  the  present  problem. 
Moreover,  it  is  easy  to  show  with  the  help  of  (2.14),  (2.10),  (2.1 1)  for  k=2,  (2.6)  and  (2.7)  that,  at 
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any  phase  boundary  with  high  strain  on  the  left,  low  strain  on  the  right,  the  driving  traction  is 
given  by 


(2.16) 


As  in  the  analyses  in  [2,3],  a  kinetic  relation  is  to  be  prescribed  at  a  phase  boundary;  we 
take  it  to  have  the  form  of  a  relation  between  driving  traction  f  and  phase  boundary  velocity  s: 

f(t)  =  (p(s(t)),  (2.17) 

where  (p  is  a  function  determined  by  material  2.  It  is  assumed  that  (p(s)  is  a  continuous  function 
that  increases  monotonically  with  s.  The  entropy  inequality  (2.15)  imposes  the  restriction 

{p(s)s>0  (2.18) 

on  the  kinetic  response  function  tp;  this  and  the  continuity  of  (p  imply  in  panicular  that  (p(0)=0. 

Again  following  the  arguments  in  [2,  3],  we  impose  a  nucleation  criterion  for  the 
initiation  of  a  phase  transformation  from  the  low-strain  phase  to  the  high-strain  phase  in  material 
2:  such  a  transformation  takes  place  through  the  emergence  at\  =  0of  a  phase  boundary 
whenever  the  associated  driving  traction  f  would  be  at  least  as  great  as  a  given  critical  value  f„ 
th^T  is  also  determined  by  the  nature  of  material  2.  After  entering  the  bar,  the  phase  boundary 
moves  to  the  right  in  accordance  with  the  kinetic  relation  (2.17).  We  assume  that  the  critical 
value  of  driving  traction  satisfies 


0  <  f  <  —  (a^ . 
*  2  M 


a  ) . 
m 


(2.19) 
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The  lower  bound  in  (2.19)  is  a  necessary  consequence  of  the  entropy  inequality  (2.15);  the  right 
inequality  guarantees  that  material  2  will  support  slowly  propagating  phase  boundaries  as  well  as 
the  fast  ones  permitted  here;  see  the  discussion  of  quasi-static  phase  transitions  in  [1-3]. 
Imposing  the  upper  bound  on  f^  in  (2.19)  also  simplifies  the  details  of  nucleation  in  the  problem 

to  be  treated  here,  so  we  adopt  it  even  though  it  is  not  strictly  necessary  to  do  so. 

Wave  propagation  properties  of  the  material  1/material  2  bar  modeled  above  are 
ultimately  to  be  compared  with  corresponding  properties  of  a  material  1/material  3  bar,  in  which 
material  3  is  a  linearly  elastic  material  whose  density  and  Young’s  modulus  coincide  with  their 
counterparts  in  the  low-strain  phase  of  material  2:  =  p^,  =  ^2*  Thus  to  treat  the  material  1/ 

material  3  bar,  one  must  modify  the  basic  field  equations  (2.1)-(2.6)  and  jump  conditions  (2.8), 

.  A  A 

(2.9)  only  to  the  extent  of  replacing  G^(y)  by  <s^{y)  in  (2.4),  and  thereupon  replacing  (2.6)  by 

^3(7)  =  ^2^.  (2.20) 

Equations  (2.10)-(2.15)  pertaining  to  the  energetics  of  the  bar  remain  valid  but  trivial,  since  for 
discontinuities  in  either  material  1  or  material  3,  (2.14)  always  yields  f  =  0  in  place  of  (2.16). 
Since  phase  transitions  cannot  occur  in  either  material  1  or  material  3,  the  kinetic  relation 
relation  (2.17),  the  attendant  restriction  (2.18)  and  the  nucleation  criterion  all  become  irrelevant 
and  are  discarded. 

3.  The  wave  propagation  problem.  We  now  formulate  the  wave  propagation  problem 
to  be  considered.  We  suppose  that,  at  time  t  =  -  <»,  an  incident  wave  bearing  a  given  tensile 
strain  Yj  >  0  and  a  given  particle  velocity  v^  is  initiated  at  x  =  -  <»  in  material  1,  traveling  to  the 

right  and  striking  the  bimaterial  interface  x  =  0  at  time  t  =  0-.  At  t  =  0+,  a  sound  wave  will  be 
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reflected  from  the  interface  back  into  material  1,  and  a  second  sound  wave  will  be  transmitted 
into  the  right  half  of  the  bar,  which  heretofore  was  at  rest  and  unstrained,  and  therefore  in  the 
low-strain  phase  of  material  2.  If  the  incident  wave  fails  to  be  strong  enough  to  nucleate  a  phase 
transition  in  material  2,  only  these  two  sound  waves  will  be  generated.  On  the  other  hand,  if 
nucleation  does  occur,  there  will  also  be  a  phase  boundary  that  emerges  from  the  interface  x  =  0 
and  travels  to  the  right,  leaving  the  particles  of  material  2  that  are  behind  it  in  the  high-strain 
phase. 


Thus  we  are  given  an  incident  wave  of  the  form 


7(x,  t)  = 


Yr 

0, 


X<Cjt, 

X>Cjt, 


V(x,t)  =  ^ 


0, 


X<Cjt, 

X>Cjt, 


for  t  <  0, 


(3.1) 


where  Yj  and  Vj  are  given  constants,  and  Cj  is  the  sound  speed  in  material  1.  Observe  that, 
because  Y  and  v  are  piecewise  constant,  the  field  equations  (2.1)  -  (2.5)  are  trivially  satisfied 
away  from  x  =  Cjt  for  t  <  0.  When  one  applies  the  jump  conditions  (2.8),  (2.9),  specialized 

appropriately  for  material  1,  to  the  discontinuity  at  x  =  Cjt  for  t  <  0  ,  one  finds  that 

Cj  =(lti/pj)^'^,  Vi  =  -C]Yj,  (3.2) 

thus  determining  the  sound  speed  in  material  1  and  imposing  a  restriction  on  the  strain  and 
particle  velocity  in  the  incident  wave. 

Given  the  incident  wave,  and  assuming  first  that  nucleation  does  occur,  we  seek  functions 
Y(x,  t)  and  v(x,  t)  of  the  following  form  on  the  upper  half  of  the  x,  t-plane  (see  Figure  2); 
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Vj,  X<-Cjt, 

''r.  -Cjt<x<0, 

,  0  <  X  <  St ,  for  t  >  0.  (3.3) 

,  St  <  X  <  C2t , 

0  ,  C2t  <  X  <  oo  , 

+  + 

here  the  two  constants  y_,  v„  associated  with  the  reflected  wave  and  the  four  constants  Y_  ,  v  of 
the  transmitted  wave  are  to  be  determined,  as  are  the  respective  constant  speeds  s  and  c^  of  the 
phase  boundary  and  the  sound  wave  in  the  low-strain  phase  of  material  2.  It  is  assumed  in  (3.3) 
that  s  is  less  than  C2,  so  that  the  phase  boundary  moves  subsonically  in  material  2;  one  can  show 
that  this  is  in  fact  necessary. 

Since  Y  and  v  of  (3.3)  are  piecewise  constant,  the  field  equations  (2.1)-(2.6)  are  trivially 
satisfied  away  from  discontinuities.  We  shall  speak  of  the  problem  of  determining  the  unknown 
constants  in  (3.3)  from  the  jump  conditions  (2.8),  (2.9)  and  the  entropy  inequality  (2.15)  at  the 
two  sound  waves,  the  bimaterial  interface  and  the  phase  bound  ry  as  the  wave  propagation 
problem.  We  shall  find  that  this  problem  has  a  1 -parameter  family  of  solutions,  parameter  s;  s  is 
then  determined  by  enforcing  the  kinetic  relation  at  the  phase  boundary  x  =  st. 

If  the  incident  wave  fails  to  nucleate  a  phase  transition  in  material  2,  the  appropriate  form 
for  the  solution  corresponding  to  the  incident  wave  (3.1)  is  that  obtained  by  deleting  from  (3.3) 
the  wedge  0  <  x  <  st  associated  with  the  transformed  material: 

00  <  X  <  -  Cjt  , 

-Cjt<x<0, 

for  t  >  0.  (3.4) 

0  <  X  <  c.,t , 

C2t  <  X  <  00  , 
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The  form  (3.4)  is  also  appropriate  for  the  composite  bar  composed  of  materials  1  and  3. 
Formally,  one  can  infer  properties  of  the  solution  whose  form  is  (3.4)  from  that  of  the  form  (3.3) 
by  setting  s  =  0  prior  to  enforcing  the  kinetic  relation,  thus  avoiding  the  need  to  consider  (3.4) 
separately  from  (3.3). 

4.  Determination  of  y  and  v. 

Case  1.  Nucleation  occurs.  We  first  assume  that  the  incident  wave  causes  a  phase 
transition  to  occur  in  material  2,  and  we  construct  the  corresponding  reflected  and  transmitted 
waves. 


When  one  imposes  the  jump  conditions  (2.8),  (2.9)  at  the  reflected  sound  wave  x  =  -  c^t 
in  material  1,  at  the  bimaterial  interface  x  =  0,  at  the  phase  boundary  x  =  st  and  at  the  transmitted 

sound  wave  x  =  C2t  in  material  2,  one  obtains  a  system  of  eight  algebraic  equations  for  the  eight 
i  i 

quantities  v^^,  y^,  v.^.,  and  s.  Of  these  equations,  the  two  arising  from  the  reflected  sound 
wave  are  not  independent  because  of  (3.2)  j,  and  the  two  arising  from  the  transmitted  sound  wave 
determine  C2  as 


=  (4.1) 


and  then  reduce  to  a  redundant  pair.  This  leaves  six  independent  equations  for  the  seven 
i  i 

unknowns  y^,  v^^,  y.^,  v.^,  and  s  .  These  equations  are  readily  solved  for  the  y’s  and  v’s  in  terms 
of  s,  furnishing 


1+|3  ^i'a(l+P)  C2  +  S  ^0’ 


(4.2) 
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2aB  s  +  (l-Kp)c 

1-^P  ^^‘(1  +  P)(C2-hs)^«’ 

(4.3) 

+  2ap  s  (pt2  +  s) 

(4.4) 

2p  1  s 

''r  ''t-  1+pVl  a(l+P)  Cj  +  s  ^1^0’ 

(4.5) 

+  2  p  (Pco 

^  1  +  P  (1 -(-P)(c5-s^) 

(4.6) 

where  we  have  introduced  the  symbols 


a  =  Cj/c2,  P  =  p,c,/p2C2, 


(4.7) 


for  the  respective  ratios  of  the  sound  speeds  and  the  mechanical  impedances  of  materials  1  and  2, 
and  we  have  also  eliminated  the  particle  velocity  Vj  of  the  incident  wave  from  the  results  by 
using  (3.2)2- 

The  representations  (4.2)-(4.6)  for  the  strains  and  particle  velocities  involve  the  as  yet 
unknown  value  of  the  phase  boundary  speed  s,  which  will  ultimately  be  determined  by  the 
kinetic  relation. 

To  assure  that  arid  are  respectively  in  the  high-  and  low-strain  phases  of  material  2, 
we  must  necessarily  enforce  the  phase  segregation  inequalities  y^.  >  y^,  y^.  ^  Yj^  .  By  (4.3)  and 
(4.4),  these  inequalities  are  equivalent  to 


1  -hP 

Yr^GJs)  =  ^^  Y 


1  (1  +  p)c  +  s 


m 


2aP  ‘tn  2ap 


C2  +  S 


^0' 


0  <  s  <  c, 


2  ’ 


(4.8) 
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l+p 


1  s  (Pc  +  s) 


Yl  <  Gj^(s)  -  ^2  _  .^2  ^0  ’ 


0  <  s  <  c, 


2  • 


(4.9) 


The  functions  and  G^  defined  in  (4.8),  (4.9)  are  both  monotonically  increasing  with  s; 
moreover,  Gj^(O)  >  0,  and  Gj^(s)  tends  to  +<»  as  s  tends  to  C2,  while  G^(C2)  is  finite.  The 
inequalities  in  (4.8),  (4.9)  are  restrictions  on  the  datum  Yj  and  the  phase  boundary  speed  s.  These 
restrictions  are  illustrated  in  the  s,  Yj-plane  of  Figure  3,  where  the  curves  Yj  =  Gj^(s)  and 
Fj^:  Yj  =  Gj^(s)  are  shown  schematically;  in  the  figure,  it  has  been  assumed  for  definiteness  that 
G^(0)  >  0  (or  equivalently  >  0),  though  this  need  not  be  the  case.  Only  the  pairs  (s,  Yj)  that 
correspond  to  points  on  or  between  these  two  curves  are  permitted  by  the  phase  segregation 
requirement. 

To  find  the  driving  traction  f  at  the  phase  boundary,  one  substitutes  for  Yj  and  from 
(4.3),  (4.4)  into  (2.16);  the  result  is 


4aP 

TTp 


r(l+P)c5-2Pc2S-2k^l 

- — -TTo - Yq-^M 


(1  +  P)(c?  -  s^) 


'm 


(4.10) 


The  curve  F^  in  the  s,  Yj-plane  along  which  f  =  0  is  therefore  given  by 


(4.11) 


we  call  Fq  the  Maxwell  curve.  For  Okskc^,  the  right  side  of  (4. 1 1)  is  a  monotonically 
increasing  function  of  s  that  tends  to  +«»  as  s  tends  to  C2.  Also,  from  (4. 1 1),  (4.8),  (4.9),  (2.6)  and 
the  assumption  that  the  Maxwell  stress  is  positive,  one  has 
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(4.12) 


1  +  P  1  P  ^0 


The  curve  is  also  shown  schematically  in  Figure  3  .  The  requirement  f  >  0,  imposed  by  the 

entropy  inequality  (2.15),  holds  only  for  points  (s,  y^)  in  the  closed  curvilinear  strip  between  F^ 

and  F- .. 

M 

Finally,  we  must  assure  that  y(x,  t)  >  -1  everywhere  in  the  x,t-plane.  One  can  show  that, 
if  (s,  Yj)  lies  in  the  strip  between  F^  and  Fj^,  the  strains  y^^  and  y^.  defined  in  (4.2),  (4.3) 

automatically  fulfill  this  requirement.  On  the  other  hand,  it  turns  out  that  y.^,  >  -1  if  and  only  if 
(s,  Yj)  lies  above  the  curve  F  ^  defined  by 


■r 


1+P 


2ap 


s  (pc,  +  s) 

— 


0  <  s  <  c„ 


C2-S 


(4.13) 


From  (4.13)  and  (4.9),  it  can  be  seen  that  F^^  always  lies  above  F  y  After  some  algebra,  one 
shows  that  F  ^  intersects  F^  exactly  once;  the  intersection  occurs  at  a  value  s  =  s^  given  by 


1/2 


C2  <  C2. 


(4.14) 


For  values  of  s  to  the  left  of  s,^,  F_j  lies  below  F^,  while  for  s„  <  s  <  C2,  the  reverse  is  true. 

Let  5  stand  for  the  closed  curvilinear  strip  bounded  above  by  F^^  and  below  in  pan  by  F^ 
and  in  part  by  F  j;  S  is  shown  hatched  in  Figure  3.  It  then  follows  that  for  each  value  of  the 
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incident  strain  that  lies  in  the  interval  (0^(0),  +<»),  there  is  a  1 -parameter  family  (parameter  s) 
of  admissible  fields  7(x,  t),  v(x,  t)  of  the  form  (4.2)-(4.6)  that  fulfill  all  of  the  conditions  of  the 
wave  propagation  problem.  For  each  given  Yj  in  this  interval,  the  permissible  range  of  the 
parameter  s  is  that  corresponding  to  the  end-points  of  the  associated  horizontal  line  segment  in 
the  s,  Yj-plane  connecting  either  or  the  vertical  axis  on  the  left  to  either  F^  or  F  j  on  the  right, 

as  appropriate.  Each  of  these  fields  that  corresponds  to  a  positive  value  of  s  involves  a  phase 
transition.  For  each  value  of  Yj  outside  (0^(0),  -h»),  there  is  no  solution  to  the  wave  propagation 

problem  that  involves  a  phase  transition;  as  we  shall  see  in  Case  2  below,  the  only  solution 
available  for  such  an  initial  datum  Yj  is  one  in  which  no  phase  transition  occurs. 


Case  2.  Nucleation  does  not  occur.  In  this  case,  we  can  find  the  solutions  of  the  form 
(3.4)  arising  from  the  incident  wave  (3.1)  by  formally  setting  s  =  0  in  the  expressions  (4.2)-(4.6). 
This  yields 


Y(x,  t)  = 


—  y  -  c  t  <  X  <  0, 

1+P  ^  ^  2ap 

v(x,  t)  =  -  - — -  c  Yt  ,  -  c  t  <  X  <  c-t ,  t  >  0  ,  (4.15) 

Zap  „  I  .HP  21  1  2 

—  y,,  0<x<Cjt, 


for  the  reflected  and  transmitted  waves  in  the  absence  of  the  phase  transformation.  Since  the 
solution  (4.15)  involves  only  sound  waves,  the  entropy  inequality  (2. 15)  holds  automatically. 
The  phase  segregation  requirement  demands  that  y(x,  t)  <  Yj^  for  0  <  x  <  c.,t;  from  (4. 15)  and 
(4.9),  this  leads  to 


(4.16) 


For  each  positive  initial  datum  Yj  satisfying  (4.16),  the  fields  (4.15)  comprise  a  solution  of  the 
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wave  propagation  problem  that  is  uniquely  determined  by  the  data  and  in  which  no  phase 
transition  takes  place.  Since  the  incident  strain  Yj  has  been  assumed  positive  (tensile)  throughout, 

the  requirement  7(x,  t)  >  -1  is  automatically  satisfied  by  the  strain  field  of  (4.15)j. 

It  is  important  to  note  that  (4.15)  also  represents  the  solution  to  the  wave  propagation 
problem  in  which  the  left  half  of  the  bar  is  composed  of  material  1 ,  while  the  right  half  is  made 
of  material  3.  In  this  case,  however,  (4.15)  represents  a  solution  for  all  tensile  values  of  the 
initial  datum  Yj,  and  not  merely  for  those  satisfying  (4. 16). 

Combining  the  conclusions  reached  in  Cases  1  and  2  establishes  the  following  results. 

(i)  For  each  value  of  the  initial  datum  Yj  in  the  interval  (0,  0^(0)),  there  is  exactly  one  solution  to 
the  wave  propagation  problem;  it  is  given  by  (4.15)  and  does  not  involve  a  phase  transition,  (ii) 
For  each  Yj  in  (0^(0),  Gj^(O)],  there  is  a  1 -parameter  family  of  solutions  (4.2)-(4.6)  to  the  wave 

propagation  problem  that  involve  a  phase  transition,  and  a  single  solution  (4.15)  that  does  not. 
(iii)  For  values  of  Yj  >  Gj^(O),  there  is  only  the  1-parameter  family  of  soutions  (4.2)-(4.6),  each 
of  which  involves  a  phase  transition. 

The  ambiguity  remaining  when  Yj  is  in  (0^(0),  Gj^(O)]  (case  (ii)  above)  is  resolved  by 
first  invoking  the  nucleation  criterion,  then  the  kinetic  relation.  Setting  f  in  (4. 10)  equal  to  the 
nucleation  value  f*  defines  a  curve  F,:  =  G^(s)  in  the  s,  Yj -plane;  we  omit  the  formula  for  G*, 

and  we  do  not  show  F^  in  Figure  3.  The  constitutive  inequality  (2.19)  can  be  shown  to  guarantee 
that  F^  always  lies  in  the  between  F^  and  Fj^,  and  that  Gp(0)  <  G^(0)  <  Gj^(O).  We  assume 
further  that  f^  is  such  that  F^  lies  in  the  slightly  smaller  strip  S  in  Figure  3.  At  points  on  or 
between  F,  and  F^^,  (4.10)  yields  values  of  f  at  least  as  great  as  f^.  It  then  follows  that 
nucleation  will  occur  if  the  incident  strain  Yj  is  at  least  as  great  as  the  critical  value  Yj  given  by 
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(4.17) 


Yi 


When  the  incident  wave  carries  a  strain  Yj  that  is  less  than  Yj  .  no  phase  transition  occurs,  and  the 

appropriate  solution  to  the  wave  propagation  problem  is  that  given  by  (4. 15).  On  the  other  hand, 

* 

if  Yj  is  at  least  as  great  as  Yj  .  nucleation  will  take  place  in  material  2,  a  phase  boundary  will 
emerge  at  x  =  0,  and  the  appropriate  solution  must  be  selected  from  the  1 -parameter  family  of 
admissible  fields  (4.2)-(4.6)  involving  a  phase  transition.  The  kinetic  relation  (2.17)  provides  the 
mechanism  for  this  selection.  Substituting  for  f  from  (4. 10)  into  (2. 17)  yields  the  equation  for  s: 

Yj  =  G  (s)  =  G  (s)  +  —  (pis) ,  y  >y  (4. 1 8) 


where  (p  is  the  kinetic  response  function  of  the  material.  Because  (p(s)  and  Gq(s)  both  increase 

monotonically  with  s,  the  same  is  true  of  the  right  side  of  (4.18).  Since  (p(0)  =  0,  one  has 
G|^(0)  =  Gq(0);  moreover,  G|^(s)  tends  to  -t-<»  as  s  tends  to  c^,  and  0  <  Gq(s)  <  Gj^(s)  <  Gj^(s)  for 
0  <s<c^.  Thus  the  curve  in  the  s,  Yj-plane  represented  by  (4.18)  lies  between  F^  and  Fj^^ . 
We  further  assume  that  tp  is  such  that  Fj^  lies  in  the  strip  S  (see  Figure  3).  Clearly  (4.18) 
determines  exactly  one  value  of  s  for  each  given  y^  in  (Gp(0),  +«).  When  inserted  into 

(4.2)-(4.6),  this  value  of  s  completes  the  determination  of  the  response  to  the  incident  wave  when 
a  phase  transition  is  nucleated. 

5.  Reflectivity.  For  a  given  incident  wave  in  material  1 ,  we  wish  to  compare  the 
strengths  of  the  reflected  waves  when  (i)  the  right  half  of  the  bar  consists  of  material  2,  and  (ii) 
the  right  half  is  composed  of  material  3.  To  motivate  our  notion  of  the  strength  of  the 
transmitted  or  reflected  wave,  it  is  helpful  to  reconsider  the  energetics  of  the  composite  bar. 
Consider  first  the  case  of  the  material  I /material  2  bar.  By  extending  the  energy  considerations 
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of  Section  2  to  a  portion  of  the  bar  long  enough  to  include  the  four  discontinuities  at  x  =  -  c^t, 

X  =  0,  X  =  St  and  x  =  C2t,  one  can  establish  the  following  identity  among  energy  rates  for  fields  of 
the  form  (3.3): 

Oj  =  +  e.p  +  fAs ,  (5.1) 


where 


e 


T 


W2(V 


{r?(YR-Yi)^^(vR-Vj)2}p^CjA, 

+  i  jp^sA  + 1 1  ^  i  I 


(5.2) 

(5.3) 

(5.4) 


f  is  the  driving  traction  at  the  phase  boundary,  A  is  the  cross-sectional  area  of  the  bar  and  s  is  the 

i  i 

phase  boundary  speed.  In  (5.3),  (5.4),  y.^,,  Vj^  and  v.^,  are  given  in  terms  of  s  by  (4.2)-(4.6), 

and  W2(yj,)  refers  to  the  stored  energy  density  of  material  2  .  We  call  Cj,  Cj^  and  e.^,  the 
incident,  reflected  and  transmitted  energy  rates;  note  that  all  have  the  units  of  energy  per  unit 
time,  and  all  are  non-negative.  By  (2.15),  the  dissipation  rate  fAs  is  non-negative  as  well.  Thus 
neither  e^^  nor  e.^,  is  greater  than  ej . 

Setting  s  =  0  in  (5.  l)-(5.4)  provides  the  corresponding  identity  among  energy  rates  for  the 
material  1/material  3  bar  in  which  no  phase  transformation  can  occur. 

The  reflectivity  and  transmissivity  are  defined  by 
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(5.5) 


respectively;  they  are  measures  of  the  strengths  of  the  reflected  and  transmitted  waves  relative  to 
the  incident  wave.  We  are  interested  in 


Qr(s,7i)  = 


Qj(^,  Yj)  = 


(5.6) 


which  are  the  ratios  of  reflectivity  and  transmissivitiy  in  the  material  1/material  2  bar  to  their 
respective  counterparts  in  the  material  1 /material  3  composite.  We  study  Q„  in  the  present 
section,  Q.^,  in  the  next.  When  p  =  1,  the  reflectivity  <lR|j_Q  of  the  material  1/material  3  bar 
vanishes,  and  (5.6)  j  fails  to  define  we  exclude  this  possibility  for  simplicity  by  assuming 
that  P  5*  1. 

The  reflectivity  ratio  is  a  function  on  the  admissible  strip  S  shown  hatched  in  Figure 
3.  For  a  given  kinetic  relation,  only  the  values  of  at  points  on  the  associated  kinetic  curve  F 

K  k 

are  relevant;  for  a  given  value  of  the  strain  y.  in  the  incident  wave,  only  the  value  of  Q„  at  the 
point  on  whose  ordinate  is  Yj  is  relevant.  If,  for  example,  one  has  Qr(s,  Yj)  <  1  this  point, 
then  the  reflectivity  of  the  material  1/material  2  bar  is  less  than  that  of  the  material  1/material  3 
bar,  so  that  reflection  has  been  diminished  by  the  occurrence  of  the  phase  transition.  We  study 
the  properties  of  as  a  function  on  5. 

From  (5.6)j,  (5.5)j,  (5.3),  (4.2)  and  (4.5),  one  finds 
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Qr(s,  Y,)  =  i 


1  -• 


1 


a(l-3)  Yj  +  k 


(5.7) 


it  may  be  noted  from  Figure  3  that  Yj  ^  >  0  at  all  points  on  5;  this  fact,  together  with  the 

assumption  that  p  ^  1,  guarantees  that  Qj^  is  well  defined  on  S.  Of  course,  Qr(0,  Yj)  =  we  now 
determine  the  locus  of  points  (s,  y.)  in  the  s,  Y,-plane  at  which  Q  =  1  and  s  0;  such  points 
correspond  to  situations  in  which  the  reflectivity  is  unaltered  by  the  phase  transition.  This  locus 
consists  of  the  interior  points  on  the  curve  defined  by 

K 


Yt  = 


1 


I  2a(l-P)  +  s  '0 


-y  0  <  s  <  c., . 


(5.8) 


Only  those  points  on  F^,  if  any,  that  lie  in  the  admissible  strip  S  are  relevant.  If  a  portion  of  F^ 
lies  in  the  interior  of  S,  the  strip  is  divided  into  two  pans:  on  one,  Qj^  <  1,  while  on  the  other, 

Q  >  1 .  Whether  F^  has  points  in  common  with  the  interior  of  S  is  determined  by  two  material 

IX  ix 

parameters  P  and  X,  where  P  =  (positive)  impedance  ratio  introduced  earlier,  and 

X  is  defined  by 


Y  +  Y 

X  =  -^ - E  .  (5.9) 

^0 

By  (2.7)  and  our  assumption  that  the  Maxwell  stress  is  positive,  we  have  A.  >  1.  We  first 
show  that  F^/ai/^  to  intersect  the  interior  of  S  if  either  of  the  following  conditions  holds: 

(i)  P  >  1  or  (ii)  0  <  P  <  1  and  X  >  A(P^),  (5.10) 
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where 


A(z)=l+-^ 

2 


1  +  z-[(l  -z)(l  -h3z)] 
1  -  z 


1/2  1 


y ,  0  <  z  <  1. 


(5.11) 


One  can  show  that  A(z)  increases  monotonically  from  the  value  1  at  z  =  0,  tending  to  +<»  as  z 
tends  to  1  . 

To  establish  (5.10),  we  first  note  that  by  (4.1 1),  (4.12)  and  (5.8),  the  curve  lies  below 
the  Maxwell  curve  near  both  endpoints  of  the  interval  0  <  s  <  c^  ,  so  that  F^  cannot  have 
points  in  common  with  the  interior  of  the  admissible  strip  S  unless  f'  intersects  F-  at  least  twice. 
To  investigate  such  intersections,  one  equates  the  right  sides  of  (4.1 1)  and  (5.8),  obtaining  a 
quadratic  equation  for  s.  It  is  readily  shown  that  (5.10)  describes  precisely  the  conditons  under 
which  this  equation  has  no  real  roots  in  the  interval  (0,  c^).  Thus  (5.10)  is  indeed  sufficient  to 

assure  that  F^  fails  to  intersect  the  interior  of  S.  Conversely,  when 

|3  <  1  and  X  <  A(P^),  (5.12) 

the  quadratic  equation  mentioned  above  has  two  real  roots  in  (0,  c^),  corresponding  to  two  points 
at  which  F^  and  F^  intersect.  By  a  calculation  too  lengthy  to  be  included  here,  one  can  show 
that  the  smaller  of  these  two  roots  lies  in  (0,  s^)  (Figure  3),  so  that  at  the  left  point  of  intersection, 
F^  indeed  enters  5. 

We  are  now  in  a  position  to  delineate  conditions  under  which  the  occurrence  of  the  phase 
transition  always  increases  the  reflectivity,  always  decreases  it  or  might  do  either.  Figure  4 
shows  three  open  regions  marked  I,  II  and  III  in  the  quadrant  P  >  0,  X  >  1  in  the  p,  plane. 
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Region  IE  corresponds  precisely  to  the  inequalities  (5.12)  and  hence  to  the  case  in  which 
enters  S.  It  follows  that  the  reflectivity  ratio  Qj^  is  less  than  one  at  some  points  (P,  X)  in  this 
region,  greater  than  one  at  others;  the  actual  effect  obtained  depends  on  the  particular  kinetic 
relation  and  possibly  on  the  panicular  value  of  the  strain  in  the  incident  wave  as  well.  The 

region  marked  I  corresponds  to  values  of  P  and  X  for  which,  on  the  interior  of  the  admissible 
strip  S  of  Figure  3,  reflectivity  at  the  interface  is  always  less  in  the  material  1/material  2  bar  than 
it  is  in  the  material  1/material  3  bar.  Thus  assuming  that  nucleation  has  occurred,  the  phase 
transformation  induced  in  material  2  by  the  incident  wave  always  acts  to  reduce  reflection  if 
(P,  X)  corresponds  to  a  point  in  region  I,  regardless  of  the  particular  kinetic  relation  involved  .  In 
region  II,  precisely  the  reverse  is  true;  reflectivity  is  always  increased  at  points  in  the  interior  of 
5  by  the  phase  transformation.  For  values  of  (P,  X)  belonging  to  the  various  regions  I,  II  and  III, 
Figure  5  schematically  shows  the  curve  F*  in  the  s,  y -plane. 

If  the  impedance  of  material  1  is  greater  than  that  of  the  parent  phase  of  material  2,  then 
P  >  1,  and  it  follows  from  (5.7)  that  Qj^(s,  y^  increases  with  s  for  each  fixed  yj .  Thus  at  each 
given  y ,  the  maximum  value  of  the  reflectivity  ratio  Q,,  occurs  on  the  right  boundary  of  5,  while 
the  minimum  value  occurs  on  the  left  boundary.  Thus  when  P  >  1 ,  fast  kinetics  promote  the 
increase  of  reflectivity,  while  slow  kinetics  are  best  for  reducing  it. 

If  the  impedance  ratio  P  <  1 ,  the  behavior  of  Q_  as  a  function  of  phase  boundary'  speed  s 
at  fixed  incident  strain  y,  is  more  complicated;  under  certain  circumstances,  may  vanish  at 

1  K 

some  points  in  S,  corresponding  through  (5.7)  to  a  minimum  in  Qj^. 
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6.  Transmissivity.  We  state  here  without  proof  some  results  pertaining  to  the 
transmissivity  ratio  Q^,  whose  behavior  is  more  complicated  than  that  of  the  reflectivity  ratio 

Using  (5.6)2,  (5.5)2,  ‘ 


8a2p2  I  Yj  ) 


T,(s), 


(6.1) 


where 


Tj(s)  =  [2p  +  (1  +  P)k/C2] 


C2  +  S 


(6.2) 


c  s 

T  (s)  =  (X  -  1)(1  +  P)2s/c2  +  -  T  (s), 

^  ^  (c^-s^)(C2  +  s)  ^ 


(6.3) 


•2 

T,(s)  =  2p2  +  (l+4p  +  p2)l  +(l-p2)L  , 


(6.4) 


and  X  is  defined  in  (5.9).  It  is  possible  to  show  that,  for  P  >  1,  one  has  Q.^,  <  1  at  all  interior 
points  of  the  admissible  strip  5;  on  the  other  hand  if  P  <  1,  one  finds  that  Q.^.  takes  values  greater 
than  1  at  some  points  in  5,  less  than  1  at  others. 
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